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ABSTRACT 


Equations are derived for finding the stiffness matrix for a twelve- 
node shear panel undergoing parabolic deformation. The deflections of a 
hypothetical cantilever beam composed of five twelve-node shear panels 
are then compared with those of six other equivalent cantilever beams. 

Two are classical strength of materials beams, one excluding shear 
effects, and the other including shear effects. The third is a theory 
of elasticity beam. The last three beams are composed of eight-node 
linear deformation shear panels, one composed of one row of five panels, 
the second composed of two rows of five panels per row, and the third 
composed of three rows of five panels per row. Tip deflection of the 
five twelve-node shear panels beam is 93 percent of the strength of 
materials beam excluding shear, whereas that of the best eight-node 
beam is 73 percent. 


NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 



Symbol 
A, B, T 


DEFINITION OF SYMBOLS 


Definition 


C e* C f’ C g > °h 
B-g, Bg j Bji 

E 

® > f j § 3 h 


length, width, and thickness of a shear panel 
interval between the S- and Y-axes 

interval between the R- and X-axes 

modulus of elasticity 

denote the left, right, lower, and upper surfaces of 
a shear panel 



flexibility matrix element defined as the deflection 
at i caused by a unit force at j 


V K f> v 


K 


h 


latus rectum of a parabolic surface 


P force at 1 

i 

R, S Cartesian coordinate system with one rotational and 

two translational degrees of freedom and origin 
attached to and tangent to the vertex of a parabolic 
surface 


R e ’ R f’ R g» 
S e ’ S £’ Sg’ 



coordinates of the end points of a parabolic surface 
in the R,S system 


R 0’ R p’ R h 

S 0’ S p> S tJt’ S h 


Coordinates of the point of intersection of a 
S 10 ~, Six"j or p - vector with a respective e-, f- , 

g-, or h -surface 


S . 


stiffness matrix element defined as the force at i 
caused by a unit deflection at j 


u 


strain energy per unit volume 


U 


total shear panel strain energy 


X, Y 


fixed Cartesian coordinate system 



Symbol 

5 

1 

e X’ e y 

0 , P, ijr, r\ 
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DEFINITION OF SYMBOLS (Continued) 

Definition 
deflection at i 

linear strain parallel to the X-, Y-axes 

angle of rotation of the R, S coordinate system 
relative to the X,Y system 

angular strain 

Poisson's ratio 


IV 



INTRODUCTION 


An eight-node shear panel typical of figure la has been in use 
for several years in computer matrix analysis of aircraft and missile 


4 8 4 8 



structures. However, this particular model is needlessly unrealistic 
in that it simply assumes that surfaces deform linearly as shown in 
figure lb. Actually, surface 1-3 deforms in some nonlinear shape 
f (x,y) which would be difficult to determine experimentally or derive 
theoretically without extensive recourse to the formal theory of elas- 
ticity. This paper proposes to elaborate on the model of figure 1 by 
adding four nodes and defining f(x,y) to be parabolic as in figure 2. 
The nodes 9, 10, 11, and 12 are assumed to be located initially at the 
mi d-points of their respective surfaces, and the distribution of strain 
(figure 2) remains linear as in figure lb. 
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Figure 2, Twelve-Node Shear Panel 


DERIVATIONS 


With Sg nonzero, locate the origin of the X-Y coordinate system at 
the 1-2 node, and attach an R-S coordinate system, to the vertex of the 
e-surface (figure 3). Thus, the equation of the e-surface is either 



( 1 ) 


or 


(y - 


e = 




K 


+ C 


e> 


(2a) 


where D e = ^ B, C e - 89, and K e = -B 2 /4Q g . A point on the e-surface 
is (e,y), and shear on the e-surface is 


7 e = tan" 1 (de/ dy) = tan" 1 [2(y - D e )/K fi ]. 


(3a) 


2 



X 



Progressively, nonzero 5 X and 63 introduce an angle of rotation 0 
such that 

0 = tan " 1 ~ > (4a) 


Figure 4 shows that the Sg-vector is no longer orthogonal to the e-sur- 
face hut remains always parallel to and midway between the Sx“ and 53 - 
vectors. Also, the vertex of the parabolic e-surface has rotated through 
the angle 0 while maintaining the point [ 1 / 2 (5x + 63 ) , 1/2 B] as the 
common point of intersection of (1) the 1-3 line, (2) the 8 g-vector, 
and (3) the R-axis . 


3 


fy" 8 3 “H 


FT 


)r * — S 1 — b») 
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Figure 4. Rotation of R-S Coordinate System 


Since all the five displacement vectors attached to a surface 
determine its total geometry, set & 2 and S 4 nonzero. Thus, by inspec- 
tion of figure 5, it may be observed that 


0 = tan "’ 1 


6 i ~ s 3 \ 
B + 84 - 8 


(4b) 


4 



5 



(5) 


= "2 (2Sg “ Si - 63) sm 0 


/B + 84 - 8 2 ' 
S e \ 2 cos 0 


R g - “ '2(259 “ 6 t_ - 83) cos 0 + R^. 


From geometrical similarity, R^ = R g S^, so that 


(26g - 5 ;l - S3) 3 sin 2 0 cos 3 0 

R — — — (8) 

2[(B + S 4 - S 2 ) 2 - (25 9 ~ Si - S3) 2 sin 2 # cos 2 #] 


h = ~ 


r (B + 64 - S 2 ) 2 - (2Sg - S x - S3) 2 sm 2 # cos 2 0 ] 
2 (25g “ S x - S3) cos 3 # 


The equation of the S-axis in the X-Y coordinate system is 

C e = i^ &1 + " R e jfcos & “ (y “ \ < B + § 2 + 84)^ tan 0 . (10a) 

The equation of the R-axis in the X-Y coordinate system is 

D e = -|(B + S 2 + 84) - R e sin 0 + (x - |(Si + S 3 )^) tan 0 . (11a) 


Similarly, for the f -surface of figure 6, 


p = tan" 


85 &7 

B + S s - 5s 


6 




( 13 ) 


2 " 5 g “ 87) sm p 


g - ( ^ - &6 

f \ 2 cos p 


Ef ■ " o (2&n “ 85 ” 87-) cos p + R, 


From geometrical similarity, Rp S£ = S 2 , so that 


— 65 — By) 3 sm 2 p cos 3 p 


R P = “ 


2[(B + 5 8 - 5 e) 2 “ ( 28 n - S 5 - S7) 2 sm 2 p cos 2 p] 


and 


Kf “ 


[(B + 5 8 - &s) 2 " ( 28 n " 85 - 67) 2 s in 2 p cos 2 p] 
2 (2Sn - 85 - 87) cos 3 p 


The equation of the S-axis m the X-Y coordinate system is 

C f = A + ^(65 + 87) - R £ /cos p - (y - |-(B + b s + S 6 )^) tan p 

The equation of the R-axis in the X-Y coordinate system is 

Df = j(B + 5 S + 5 6 ) - R { sin p + (x - 77(85 + 87) - tan p 


For the g -surface of figure 7 , 


= tan 


86 ~ 8g 

A + 8 S - Si 


( 14 ) 

( 15 ) 

( 16 ) 

( 17 a) 

. ( 18 a) 
( 19 a) 
(20a) 
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3 




( 21 ) 


R^ — 2 ( 25 io “ 62 “ &e) sm t); 


A +_ S5 - 5 ; 


l 2 cos ijr 


Sg — - 2 (2&io — S2 ” 5g) cos ijr 4- S 


r 


2 2 

From geometrical similarity, R g ^S g R^, so that 




(2&io “ 5 ^ - 6s) 3 sm 2 -i|r cos 3 ij; 


2[ (A + 6 5 - 5 X ) 2 - (28 10 - 8 2 - 8 6 ) 2 smSjr cos 2 ijr] 


and 


K g = - 


[(Ad- 85 “ 5 i) 2 “ (25x0 “ S 2 “ 5 6 ) 2 sin 2 ty cos 2 i|r] 
2 ( 28^0 “ &2 “ 85) cos 3 ^r 


The equation of the S-axis m the X-Y coordinate system is 


C g = “(A + 5 X + 8 5 ) - S g sin \|r - fy - ^-(5 2 + 5 6 ) ) tan ijr. 


The equation of the R.-axis m the X-Y coordinate system xs 


Dg = ■^‘(&2 + 5 e) " s g / cos t - ■|'(A + 8 X + S 5 )^ tan i|r 


Finally, for the h -surface of figure 8, 


T) - tan - 


-i („ &8 “ 

\A + 8-7 - 83 


( 22 ) 

(23) 

(24) 
(25a) 

(26a) 

(27a) 

(28a) 
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(29) 


R T) ” 2^ 2612 " 54 " s:Ln ^ 

R , = A + 87 r. s jA 

Kh \ 2 cos n / 

s h = - J (28x2 - 84 - 8g) cos T| + S^. 

From geometrical similarity, R 2 , so that 

( 25 x 2 “ 64 - 8 S ) 3 sin 2 r] cos 3 t) , 

g — ^ ... 

71 2 [(A + 87 - 83) 2 - ( 28 x 2 - 84 - & 8 ) 2 sin 2 r| cos 2 t]] 

and 

[(A + 87 - 83) 2 - (25x2 “ 5 a - S3) 2 sin 2 !] cos 2 !]] 

Kh • 

2 ( 28 j _2 “ ~ Sq) cos 3 t) 

The equation of the S-axis m the X-Y coordinate system is 

c h = ^-(A + S 3 + 8 t ) - S h sm f] - - 77(84 + 63 ) - B^) tan T|. 

The equation o£ the R-axis m the X-Y coordinate system is 

0 h = B + |(6 4 + 5 S ) " s h/ C0S R + (x - (| A + 5 3 + 5 7 )^ tan 


(30) 

(31) 

(32) 
(33a) 
(34a) 

(35a) 
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If it is reasonably assumed that deflections wrll generally be much 
less than A or B, «A or B, then 


sxn a, = ex 
cos a, = 1 
tan a = a 

s 2 = 0 

x 


s in 2 cc = 0 
cos 2 a = 1 
cos 3 a = l 


where (X = 0, p, i|r, or t], and 


Si ~ 

B + 84 - 8g 


S 5 " Sy 
B + 5 S - 8 6 


I _ S 6 - S,g 

V A + S 5 - 

5 q - 84 

^ A + Sy - S3 


(4c) 

(12b) 

(20b) 

(28b) 


K 


e 


- (B + 5 4 - S 2 ) 2 - 
2(25$ - 5i - 5 3 ) 


(9b) 


K 


f 


K 


8 


f (B + &8 - Ss) 2 -] 
_2(25 11 - o 5 - 8-7) _ 

- (A + 8 5 - Bj.) 2 - 

2(2& 10 - ~ &e)_ 


(17b) 


(25b) 
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(33b) 


(A + 5 7 - 53 ) 2 - 

2 (2 $i2 - Sa - S 8 )_ 


Cg = — (2&9 + - 5 3 ) - y0 

Cf ~ "f "2 (25i3_ + 65 - S 7 ) - yp 
c g = \ (A + &i + 5 S ) 

C h = | CA + 5 3 + 5 7 ) 

D e = \ (B + 62 + 5 a ) 

D f = | (B + S 6 + S s ) 


(10b) 

(18b) 

(26b) 

(34b) 

(Hb) 

(19b) 


D g = \ (25 10 + 62 “ 5 6 ) + xt (27b) 

D h = B + \ (2^12 + &4 - 6 S ) + XT) . (35b) 

The equatxons for each surface and its associated shear become 


(y - E e ) 2 


K +C e 

e 

(2b) 

(y - D f ) 2 


K, +G f 

(36) 


f 
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( 37 ) 


(x - C ) 2 
g = -■ k — ~ + Dg 

g 


h — 


?e = 


7f = 




C h^ 

Kh 

2(y - D e ) 

K 


+ D 

h 


- 0 


e 

2(y - D f ) 


K 

f 

2 (x - C g ) 

H 

2(x - C h ) 


P 


+ * 


+ T|. 


(38) 
(3b) 

(39) 

(40) 

(41) 


Recalling that 5^ « A or B and strain is assumed 
distributed, shear panel strains (figure 9) are 


to be linearly 


e 

x 


1 p(y - D f ) 2 

A ’ 


(y - 

K 


V® 


c f : 



(42) 


e 

y 




(x “ c g ) 


K. 


+ D 


- D. 


g 


B 


(43) 



(44) 
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Figure 9. Completely Deformed Twelve-Node Shear Panel 


The stiffness matrix elements may now be calculated by means of 
equation (45), 



16 


i 

To calculate an element, examples 1 and 2, first set all 5^, K ^ i or j, 
equal to zero in equations (42), (43), and (44); second, take the partials 
of equations (42), (43), and (44) with respect to 8^, and third, set 
8^ - 0 and 5j = 1 in the expressions for e x , Gy, y, Se x /d5^, 8ey/d8^ 5 
and 0y/c$b x . J Substitute these expressions into equation (45) and carry 
out the indicated integration. 


Example 1, Sg } x 


0 = 5i/B 


p = ijf = q - o 


K _ . . . -B 3 

e 2 (2 89 - 5i) 


K f = K g - K h = - 


C = 
e 





B 


C £ — A Dg — 0 

Cg = C h = 2 k B h = B 


e 

x 


l p2(2S 9 - &x) (y - I B ) 2 
A P 


89 





2(y - \ B) 2 

p 



Z 

B 




1 


de 


e = 5—^ = 0 

y ds 9 
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Example 2, S 1D , 9 


0 = p=Tjr=T] = 0 


K 


e 


Bf_ 

4Sg 


K 


g 


A 2 

4S 10 


K f = K h = -« 


C = 
e 


& 9 



c 


£ 


= A 


D g = S 10 


c g - Ch - 2 A D h “ B 


e 

x 


1 

A 


-45 9 (y - JB) 2 
B 2 



x 


r S 9 = 1 
&10“ 0 


I 

A 


p4(y - - B) 2 

p 


1 


x 

S8xo 


e 

y 


1 p4& l0 (X - \ A ) 2 
B A 2 




Sg “ 1 

&io = 0 


0 
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( 47 ) 



Thus, column 1, which is typical of loads caused by a corner deflec- 
tion (figure 2b), is 



12 B/A 4- 35 A/B (1 - n) 
90 


(48) 



-ET T 24 B/A - 35 A/B (1 - u) 
S 5,l - ]— 2 L Ti0 



(49) 

(50) 

(51) 

(52) 

(53) 




(54) 

(55) 

(46) 

(57) 
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(57) 


and column 9 
(figure 2a) , 


s ll, 1 


-ET 3 B/A + 10 A/E (1 - u) 
1 - [i 2 45 


S 12 } 1 


. “ E T , A+jA 
1 - p 2 V I 8 ) 


, which is typical of loads caused by a midpoint 
Is 


S l } 9 



3 B/A - 20 A/B (1 - u) 
45 



J 3,9 


ET I" 3 B/A - 20 A/B (1 - u) 

1 - |i 2 45 


S 4j 9 




^5, 9 


-ET 3 B/A + 10 A/B (1 - u) 
1 - ji 2 45 


Ss 9 

9 


ET A f iA 

1 - p 2 V 18 J 


S 7 , 9 



3 B/A 4- 10 A/B (1 - u) 
45 




(58) 

deflection 

(59) 

(60) 

1 

(61) 

(62) 

(63) 

(64) 

(65) 
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1 + ft 
18 


( 66 ) 


S 8j 9 ” 


-ET 
1 - 




S 9, 9 


ET f 24 B/A + 40 A/B (1 - u) 

1 - [a 2 [ 45 


(67) 


SlO^ 9 


ET (2 + 2(j\ 

1 “ H 2 \ 9 ) 


(47) 


S ll, 9 


-ET [~ 24 B/A - 20 A/B (1 - u) 
1 - p 2 45 


( 68 ) 


S 12, 9 


-ET / 2 + 2ia \ 

1 - ix 2 V 9 ) 


(69) 


By consideration of the geometry of the shear panel, the remaining 
120 stiffness matrix elements may be substituted in terms of those 
already calculated Figure 10 graphically illustrates these substitu- 
tions. Those with a bar over them essentially represent a rotation of 
90 degrees of the shear panel. For example. 


S 2, 2 . - S l, l/90‘ 


ET 12 A/B + 35 B/A (1 - u) 
1 - 90 


• (70) 


Actually, since the stiffness matrix is positive definite, i.e.. 


12 12 

II 


i=lj=l 


S. 


B. 5. > 0 
i 1 


= 0 


for any arbitrary 8^ 
when not all 5^ = 0, 
for all 8^=0, 


23 



j 

= 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

= 1 


2,1 

3,1 

- 6,1 

5,1 

- 4,1 

7,1 

8,1 

1,9 

- 4,9 

7,9 

- 6,9 

2 

2,1 

1,1 

- 4,1 

5,1 

- 6,1 

3,1 

8,1 

7 yi 

2,9 

3,9 

8,9 

5,9 

3 

3,1 

6,1 

1,1 

- 2,1 

7,1 

- 8,1 

5,1 

4,1 

3,9 

- 8,9 

5,9 

- 2,9 

4 

4,1 

5,1 

- 2,1 

1,1 

- 8,1 

7,1 

6,1 

3,1 

4,9 

7,9 

6,9 

— 

1,9 

5 

5,1 

4,1 

7,1 

- 8,1 

1,1 

- 2,1 

3,1 

6,1 

5,9 

- 2,9 

3,9 

- 8,9 

6 

6,1 

3,1 

- 8,1 

7,1 

- 2,1 

1,1 

4,1 

5,1 

6,9 

1,9 

4,9 

7,9 

7 

7,1 

8,1 

5,1 

- 4,1 

3,1 

- 6,1 

1,1 

2,1 

7,9 

- 6,9 

1,9 

- 4,9 

8 

8,1 

7,1 

- 6,1 

3*1 

- 4,1 

5,1 

2,1 

1,1 

8,9 

5^9 

2,9 

3,9 

9 

9,1 

10,1 

9,1 

- 10,1 

11,1 

- 12,1 

11,1 

12,1 

9,9 

- 12,9 

11,9 

- 10,9 

1 0 

10,1 

9,1 

- 12,1 

11,1 

- 10,1 

9,1 

12,1 

11,1 

10,9 

9,9 

12,9 

11,9 

1 1 

11,1 

12,1 

11,1 

- 12,1 

9,1 

- 10,1 

9,1 

10,1 

11,9 

- 10,9 

9,9 

- 12,9 

1 2 

12,1 

11,1 

- 10,1 

9,1 

- 12,1 

11,1 

10,1 

9,1 

12,9 

11,9 

10,9 

9,9 


Figure 10 Stiffness Matrix Substitutions for Twelve-Node Shear Panel 


only column 1 and the last four elements of column 9 are necessary to 
arrive at a lower triangular matrix, i S J, (although equations (59) 
to (66) serve as accuracy checks) with the upper triangular matrix given 
by Maxwell's reciprocity relationship, S = The matrix is also 

inherently singular since, to satisfy static equilibrium 

12 

) S 5=0 
i=l 


for any arbitrary j, 1 g j ^ 12, 6j f 0 
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RESULTS AND CONCLUSIONS 


Suppose, now, that the twelve-node stiffness matrix has been 
derived, and that seven cantilever beams are machined from aluminum 
with the dimensions and properties of table I 


A 

Beam (cm.) 

SOM w/o shear 10 

SOM with shear 10 

TOE 10 

5/8 2 

10/8 2 

15/8 2 

5/12 2 


TABLE I 

B T E I 

(cm) (cm) (dy/cm 2 ) (cm 4 ) u 

2 1 1/2 6.895 x 10 11 1 .3 

2 

2 

2 


1 

2/3 

2 


The first too beams (figure 11c) are analyzed using common strength of 
materials equations, excluding shear effects (SOMEXS), and including 
shear (SOMINS) . The third (also figure 11c) is analyzed from theory of 
elasticity (TOE). The next three are composed entirely of eight-node 
shear panels analyzed by matrix algebra (appendix C) , one composed of 
the first five panels of figure 11a (this beam hereafter noted by "5/8"), 
another composed of the first ten panels (10/8), and the third composed 
of all fifteen eight-node panels (15/8). The last beam is composed of 
five toelve-node panels (5/12) of figure lib. All beams are anchored 
to zero-deflection "foundations" at the left end (Sx = 67 = 833 = &x9 = 
&ss = 831 — 837 = 843 = 0 for the 5/8, 10/8, and 15/8 beams or &x - 87 = 
8x3 — &19 = &30 for the 5/12 beam). All beams are end-loaded with 
2 x 10 7 dynes of downward acting force (P^q an( ^ -^29 tor figure 11a and 
lib, respectively). Results are presented m tables II and III and fig- 
ures 12 to 16. 
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TABLE II 


SOM And TOE Cantilever Beam Deflections 


X 8 - 

(cm) 

SQM w/o Shear 
(mm) 

5 - SOM with Shear 
(mm) 

5 - TOE 
(mm) 

0 

.0000 

.0000 

.0000 

1 

-.0014 

-.0016 

-.0018 

2 

- 0054 

-.0059 

-.0062 

3 

-.0117 

-.0125 

- 0129 

4 

-.0201 

-.0210 

- 0216 

5 

-.0302 

- . 0314 

-.0321 

6 

-.0418 

- . 0432 

- . 0440 

7 

- . 0545 

- 0562 

-.0571 

8 

-.0681 

-.0700 

- 0711 

9 

-.0822 

- . 0844 

- 0856 

10 

-.0967 

-.0992 

- 1005 

Tip Rotation 

- 00145 

-.00147 

- 00149 


(Radians) 
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Figure 12. SOMEX, SOMINSj and TOE Beam Deflections 
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TABLE III 


Shear Panel Cantilever Beams External Loads and Deflections 



5/8 

10/8 

15/8 

5/12 


P. 

5. 

— t : 

5- 

P 

8. 

P. 

5- 


X 

X 

X 

X 

X 

X 

X 

X 


(dynes) 

(mm) 

(dynes) 

(mm) 

(dynes) 

(mm) 

(dynes) 

(mm) 

1 

-1 x 10 s 

.0000 

-1 x 10 s 

.0000 

-8.58 x 10 7 

.0000 

-1 x 10 s 

.0000 

2 

0 

0035 

0 

.0036 

0 

.0037 

0 

.0048 

3 

0 

0063 

0 

0065 

0 

.0066 

0 

.0084 

4 

0 

.0082 

0 

0086 

0 

0087 

0 

.0111 

5 

0 

.0094 

0 

0098 

0 

.0099 

0 

0127 

6 

0 

0098 

0 

.0103 

0 

.0104 

0 

0133 

7 

1 x 10 s 

.0000 

3.15 x 10 3 

.0000 

-4 27 x 10 T 

.0000 

-2 30 x 10 3 

0000 

8 

0 

- 0035 

0 

.0000 

0 

.0012 

0 

-.0000 

9 

0 

-.0063 

0 

- 0000 

0 

0022 

0 

.0000 

10 

0 

-.0082 

0 

0000 

0 

0029 

0 

- 0000 

11 

0 

-.0094 

0 

- 0000 

0 

0033 

0 

.0000 

12 

0 

-.0097 

0 

0000 

0 

0034 

0 

-.0000 

13 

— ■ 

- 

CO 

o 

X 

r-t 

.0000 

4.27 x 10 7 

0000 

1 x 10 s 

0000 

14 

- 


0 

-.0036 

0 

- 0012 

0 

- 0048 

15 

- 


0 

-.0065 

0 

- 0022 

0 

- 0084 

16 

- 

- 

0 

- 0086 

0 

- 0029 

0 

- 0111 

17 

- 

- 

0 

-.0098 

0 

- 0033 

0 

- 0127 

18 

w 


0 

-.0102 

0 

- 0034 

0 
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Figure 13. Eight-Node Shear Panel Beam Deflections 


30 





5 


10 

i 



Figure 14. Improvement of Eight-Node Shear Panel 


15 No. of Panels 

j . — 


O 5/8 
□ 10/8 
A 15/8 

0 5/12 for Comparison 



Maximum Deflection 







Four conclusions are drawn: (1) Of primary interest - the twelve- 

node shear panel is superior to the eight-node panel as the 5/12 beam 
tip deflection is 93 percent of the tip deflection of the SOMEXS beam, 
whereas that of the best eight-node panel beam is 73 percent (figure 15); 
(2) all the shear panel beams (figures 13 and 15) are stiffer than the 
SOM and TOE beams (figure 12); (3) increasing the number of shear panels 
composing a structure, i.e., decreasing "grain size," generally increases 
the accuracy of the results' (figure 13), but the additional accuracy 
gained by adding more rows rapidly approaches a tradeoff with the size 
of the stiffness matrix as illustrated in figure 14 and table III; and 
(4) plane sections initially orthogonal to the neutral axis tend essen- 
tially to remain plane after loading (figure 16) . 

Table IV lists the shear panel stiffness matrices in units of 
dynes/ centimeter. 
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Figure 16. Tip Rotation of 15/8 and 5/12 Shear Panel Beams 
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TA3LE 

IV 

SHEAR PANEL i 

ST1 FFilESS MATRICES 




5/8 


IB/8 


15/8 


5/12 


S[ 

1, 1]= 

.5*1 1442 3 IE 

12 

.454S1533E 

12 

.52407 051E 

12 

.46092949E 

12 

S[ 

1, 2] = 

. 18468750E 

12 

.18468750E 

12 

.184687 5 BE 

12 

.28571314E 

12 

S[ 

1, 3] = 

.56826923E 

11 

-.17248877E 

12 

-.334647 44E 

12 

•63141326E 

10 

S[ 

1, 4] = 

'.142E6731E 

11 

.14206731E 

11 

.14286731E 

11 

.45777244E 

11 

S[ 

1, 5] = 

-.31254808E 

12 

-.56826923E 

11 

.72612179E 

11 

.3157 85 1 3E 

10 

S[ 

1, 6] = 

-.14206731E 

11 

14206731E 

11 

-.142 0673 IE 

11 

-.45777244E 

11 

S[ 

1, 7] = 

-.25572115E 

12 

-.2273B769E 

12 

-.26203526E 

12 

.59983974E 

11 

S[ 

1, 8] = 

-.1846875BE 

12 

18468750E 

12 

-.1846875BE 

12 

-.2B52B833E 

11 

sc 

1, 9] = 







-.27782051E 

12 

S[ 

1,10] = 







183 10897E 

12 

S[ 

1,11]= 







-.25256410E 

12 

S[ 

1,12]= 







-.82083333E 

11 

S[ 

2, 1] = 

.18468750E 

12 

.18468750E 

12 

. 18468750E 

12 

• 2857 131 4E 

12 

S[ 

2, 2] = 

.511 44231E 

12 

.82399B38E 

12 

.11807372E 

13 

.46C92949E 

12 

S[ 

2, 3] = 

- . 14286731E 

11 

-.14206731E 

11 

-.14206731E 

11 

-.45777244E 

11 

S[ 

2, 4] = 

-.31254808E 

12 

-.72454327E 

12 

— « 1 1 1 4439 IE 

13 

.31578513E 

10 

sc 

2, 5 ] = 

. 14206731E 

11 

.14206731E 

11 

.142C6731E 

11 

. 45777244E 

11 

S[ 

2, 6] = 

.56826923E 

11 

.31254808E 

12 

.52407051E 

12 

.63141Z26E 

10 

S[ 

2, 7] = 

-.18468750E 

12 

-.1846875BE 

12 

- . 18468750E 

12 

-.20520833E 

11 

S[ 

2, 8] = 

-.25572115E 

12 

-.411995 19E 

12 

- .590 36S59E 

12 

•59983974E 

11 

S[ 

2, 9] = 







-.18310897E 

12 

S[ 

2,18]= 







-.27732051E 

12 

S[ 

2,11]= 







-.82083333E 

11 

sc 

2,12]= 







-.2525641BE 

12 

sc 

3, 1] = 

.568 269 23E 

11 

-. 17048077E 

12 

-.33464744E 

12 

.63141026E 

13 

sc 

3, 2]= 

-.14266731E 

11 

-. 14286731E 

11 

-.14226731E 

11 

-.45777244E 

11 

sc 

3, 3] = 

.51144231E 

12 

.45461538E 

12 

.52407C51E 

12 

.46092949E 

12 

sc 

3, 4] = 

-.1846875 BE 

12 

-.1846875BE 

12 

-. 1846875KE 

12 

-.28571314E 

12 

sc 

3, 5] = 

-.25572115E 

12 

-.2273C769E 

12 

-.262G3526E 

12 

.59983974E 

11 

sc 

3, 6]= 

. 18468750E 

12 

. 1846875BE 

12 

.18468750E 

12 

.2B520833E 

11 

sc 

3, 7] = 

-.31254808E 

12 

-.56826923E 

11 

♦72612179E 

11 

.3157C513E 

10 

'SC 

3, 8] = 

.142B6731E 

11 

. 14206731E 

11 

.14226731E 

11 

.45777244E 

11 

sc 

3, 9] = 







-.2778205 IE 

12 

sc 

3,10]= 







.82083333E 

11 

sc 

3,11]= 







-.25256410E 

12 

3[ 

3,12]= 







, 1831 D897E 

12 

sc 

4, 11 = 

.142 067 3 IE 

11 

.14206731E 

11 

• 1 42067 3 IE 

11 

.45777244E 

11 

sc 

4, 2] = 

-.31254308E 

12 

-.72454327E 

12 

-.11144391E 

13 

.315705 13E 

10 

sc 

4, 3] = 

- . 18468750E 

12 

-.1S468753E 

12 

-.18468750E 

12 

-.23571314E 

12 

sc 

4, 4] = 

.51144231E 

12 

.82399B38E 

12 

. 1 18G7372E 

13 

.46092949E 

12 

sc 

4, 5] = 

.18468750E 

12 

.18463752E 

12 

.1846875BE 

12 

.2B52B833E 

11 

sc 

4, 6] = 

-.25572115E 

12 

—.41 1995 1 9E 

12 

-.59036859E 

12 

.59983974E 

11 

sc 

4, 7] = 

142C6731E 

11 

-.14286731E 

11 

-»1^226731E 

11 

-.45777244E 

11 

sc 

4, S] = 

.56826923E 

11 

.31254838E 

12 

.52407051E 

12 

.63141826E 

10 

sc 

4, 9] = 







• 13318897E 

12 

sc 

4,10]= 







-.2525641 BE 

12 

sc 

4,11]= 







.82083333E 

11 

sc 

4,12]= 







-.27732051E 

12 
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TABLE IV 


[CO''T I VUED ] 




5/8 


10/8 


15/8 


5/12 


SC 

5, 1] = 

-.31254808E 

12 

-.56326923E 

11 

.72612 179E 

11 

.31 57 85 1 3E 

10 

sc 

5, 2] = 

.14206731E 

11 

.14206731E 

11 

. 142S6731E 

11 

.45777244E 

11 

sc 

5, 3] = 

-.25572115E 

12 

-.22730769E 

12 

-.26203526E 

12 

. 5998397 4E 

11 

sc 

5, 43 = 

. 184 6875 0E 

12 

.18468750E 

12 

, 18468750E 

12 

.20520833E 

11 

sc 

5, 53 = 

.51 144231E 

12 

.45461538E 

12 

.5240705 IE 

12 

.46092949E 

12 

sc 

5, 63 = 

184 6875 0E 

12 

-.18468750E 

12 

18468750E 

12 

- .2857131 4E 

12 

sc 

5, 73 = 

.56826923E 

11 

-.17&48077E 

12 

-.33464744E 

12 

.63141626E 

10 

sc 

5, 83 = 

-.14206731E 

11 

-.14206731E 

11 

14206731E 

11 

-.45777244E 

11 

sc 

5, 93 = 







-.25256410E 

12 

sc 

5,103“ 







. 18310897E 

12 

sc 

5,113 = 







-.27782051E 

12 

sc 

5,123 = 







.82383333E 

11 

sc 

6, 13 = 

-.14206731E 

11 

-.14206731E 

11 

- . 142B6731E 

11 

-.45777244E 

11 

sc 

6, 23 = 

.56826923E 

11 

•31254808E 

12 

.52407051E 

12 

.63141C26E 

10 

sc 

6, 33 = 

. 1846875 0E 

12 

.18468750E 

12 

. 18468 750E 

12 

.2G520833E 

11 

sc 

6, 43 = 

-.25572115E 

12 

- .4119951 9E 

12 

-.59036859E 

12 

.59983974E 

11 

sc 

6, 53 = 

1846875 0E 

12 

-.18468750E 

12 

-.18468750E 

12 

-.28571314E 

12 

sc 

6, 63 = 

.51 144231E 

12 

.82399038E 

12 

.118B7372F 

13 

.46092949E 

12 

sr 

6, 73 = 

. 1 420 673 IE 

11 

.14206731E 

11 

• 1 420 673 IE 

11 

.45777244E 

11 

sc 

6, 83 = 

-.31254868E 

12 

-.72454327E 

12 

-.11144391E 

13 

.3157C513E 

10 

sc 

6, 93 = 







.82C83333E 

11 

sc 

6,103= 







-.27782051E 

12 

sc 

6,113= 







. 18316897E 

12 

sc 

6,123= 







-.25256416E 

12 

sc 

7, 13 = 

- .25572115E 

12 

-.22730769E 

12 

-.26203526E 

12 

.59983974E 

11 

sc 

7, 23 = 

18468750E 

12 

-.18468750E 

12 

-.18468750E 

12 

-.20520833E 

11 

St 

7, 33 = 

-.31254808E 

12 

-. 5682692 3E 

11 

.72612 179E 

11 

. 315705 13E 

10 

SC 

7, 43 = 

-.14206731E 

11 

-.14206731E 

11 

-.14206731E 

11 

-.45777244E 

11 

S[ 

7, 53 = 

.5 6826923E 

11 

-.17048P77E 

12 

-.33464744E 

12 

.63141026E 

10 

sc 

7, 63 = 

. 1 42067 3 IE 

11 

. 142B6731E 

11 

. 1 420 673 IE 

11 

.45777244E 

11 

sc 

7, 73 = 

.51144231E 

12 

.45461538E 

12 

.5240705 IE 

12 

.46692949E 

12 

sc 

7, 83 = 

. 184687 5 0E 

12 

.18468750E 

12 

.18468750E 

12 

.28571314E 

12 

sc 

7, 93 = 







-.25256410E 

12 

sc 

7,103= 







-.82083333E 

11 

sc 

7,113= 







-. 27782051E 

12 

sc 

7,123= 







- . 1831 6897E 

12 

sc 

8, 13 = 

1846875 0E 

12 

18468750E 

12 

-.18468750E 

12 

-.26520833E 

11 

sc 

8, 23 = 

-.25572115E 

12 

-.41 199519E 

12 

-.59036859E 

12 

.599 8397 4E 

11 

sc 

8, 31= 

.14206731E 

11 

, 14206731E 

11 

• 1 426673 IE 

11 

•45777244E 

11 

sc 

8, 43 = 

.56826923E 

11 

.31254808E 

12 

.5240765 IE 

12 

.63141026E 

10 

sc 

8, 53 = 

-.14206731E 

11 

-.142C6731E 

11 

— . 1 420673 IE 

11 

-.45777244E 

11 

sc 

8, 63 = 

-.31254808E 

12 

-.72454327E 

12 

-.11144391E 

13 

.31570513E 

10 

sc 

8, 73 = 

.1 846875 0E 

12 

.18468750E 

12 

.18468750E 

12 

.28571314F 

12 

sc 

8, 83 = 

.51 144231E 

12 

.82399038E 

12 

.11807372E 

13 

.46092949E 

12 

St 

8, 93 = 







-.82083333E 

11 

St 

8,10]= 







- .25256410E 

12 

SC 

8,113= 







-.18310897E 

12 

sc 

8,123= 







-.277S2051E 

12 
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TABLE IV [CONTINUED] 


5/8 10/8 15/8 5/12 


S[ 9, 1]= 

-.277S2051E 12 

Sr 9, 2]= 

-.18310897E 12 

S[ 9, 3] = 

-.27732051E 12 

S[ 9, 43 = 

. 18318897E '12 

S[ 9, 53 = 

-.25256410E 12 

S[ 9, 63= 

.32C83333E 11 

S[ 9, 73 = 

-.25256410E 12 

St 9, 8] = 

-.82083333E 11 

S[ 9, 9]= 

.13133333E 13 

St 9,103= 

.25256410E 12 

S[ 9,113= 

-.25256410E 12 

St 9,123= 

-.25256410E 12 

SC 10, 13 = 

-.18310897E 12 

SC 1 0, 23 = 

-.27782051E 12 

S[ 1-0, 33 = 

o u2083333E 11 

S[10, 4] = 

-.25256410E 12 

SC 10, 53 = 

.18310897E 12 

St 10, 63 = 

-.27782051E 12 

St 10, 73 = 

-.82083333E 11 

SC 10, 83 = 

-.2525641SE 12 

SC 10, 93 = 

.25256410E 12 

S[ 10,103= 

. 13133333E 13 

S[ 10,113= 

-.25256410E 12 

St 10,12] = 

-.25256410E 12 

SC 11 , 13 = 

-.25256410E 12 

sen, 23 = 

-.82083333E 11 

S[11, 33 = 

-.25256410E 12 

SMI, 43 = 

.82083333E 11 

S[ 11, 53 = 

-.27782051E 12 

SC 11, 6] = 

•18318897E 12 

SC 11, 73 = 

-.27782051E 12 

SC 1 1, 83 = 

-.18310897E 12 

sell, 9]= 

-.25256410E 12 

St 11,103 = 

-.25256410E 12 

Sell, 113= 

•13133333E 13 

S[ 11,123= 

.25256410E 12 

St 12, 13 = 

-.82083333E 11 

St 12, 23 = 

“.25256410E 12 

St 12, 33 = 

.18310897E 12 

St 12, 43 = 

-.27782051E 12 

St 12, 53 = 

.82083333E 11 

St 12, 63 = 

-.25256410E 12 

St 12, 73 = 

-.18310897E 12 

Sf 12, 83 = 

-.27782051E 12 

St 12, 93 = 

-.25256410E 12 

SE 12,103= 

-.25256410E 12 

St 12,113 = 

.25256410E 12 

SC 12, 121 = 

. 13133333E 13 


*STOP* 
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APPENDIX A 

Eight-Node Shear Panel Stiffness Matrix 


By placing the origin of the X-Y coordinate system at the 1-2 node 
of the shear panel of figure la of the text, retaining the stated assump- 
tions, and setting all 8^ # 0, the strains may be easily formulated. 


Figure 17. Completely Deformed Eight-Node Shear Panel 


Thus, from figure 17, 





wher e 


s 7 ~ 


65 ~ 8^ 

is the strain at the lower surface and — — - — — at the 

A A 

upper surface. Similarly, 


e = 

y 




( 72 ) 


= ^3 - Si 5g 2 + (^ 7 " S3 - 5i \ x , /Sg - 84 5g - 5g \ ^ 

B A VB “ B ) A V A " A / B 

( 73 ) 


where 


5 3 ~ S x §6-62 


B 


+ : is the shear strain at the 1-2 node, 

A ’ 


67 - 85 S3 - 87 

^ — is the linear strain variation with x, and 


Sg — 84 Sg - S; 


is the linear variation with y. 


Again to calculate a stiffness element, S L j, set all contained 
in the expressions for e x , e y , and 7 equal to zero, K ^ r or j. Next, 
take the partials of € x , e y , and y with respect to 5.^ then set 5^ = 0 
and 5 j = 1 in the resulting expressions of e x , € y , 7, Se x /^6 1 , 5e y /^S 1} 
and hy/bb ± , and finally substitute these expressions into equation (45) 
and integrate. 

The first column of the stiffness matrix is 


c = BT 

Sl > 1 1 - p 2 


2 B/A + A/B (1 - u) 
6 


c - EI „ (1 ± U 

s *. 1 ' rri? ( s 


ET 


>3> 1 1 - p 2 


B/A - A/B (1 - u) 
'6 


( 74 ) 


( 75 ) 


( 76 ) 


40 



(77) 


S 4 i - 


s 


5,1 ~ 


> 6,1 


s 


7 , 1 ~ 


$ 8 , 1 


ET 

(1 " 3 A 

1 - y? 

^ 8 ; 

-ET 

4 B/A - A/B (1 - u) _ 

"k 

1 

T— 1 

12 

-ET 

4 - sA 

1 - f 

^ 8 J 

-ET 

”2 B/A + A/B Cl - u)~ 

1 - h 2 

12 

-ET 

(i + (A 

1 2 
1 - \r 

V 8 J • 


(78) 

(79) 

(80) 
(81) 


The geometry of substitution is identical to that of the 8x8 sub- 
matrix of figure 10 of the text 

J = 1 2 3 4 5 6 7 8 

r = 1 

2 

3 

4 

5 

6 

7 

8 


1,1 

2,1 

3,1 

- 6,1 

5,1 

- 4,1 

7,1 

8,1 

2,1 

1,1 

- 4,1 

5,1 

- 6,1 

3,1 

8,1 

7,1 

3,1 

6,1 

1,1 

- 2,1 

7,1 

- 8,1 

5,1 

4,1 

4,1 

5,1 

- 2,1 

1,1 

- 8,1 

7,1 

6,1 

3,1 

5,1 

4,1 

7,1 

- 8,1 

1,1 

- 2,1 

3,1 

6,1 

6,1 

3,1 

- 8,1 

7,1 

- 2,1 

1,1 

4,1 

5,1 

7,1 

8,1 

5,1 

- 4,1 

3,1 

- 6,1 

1,1 

2,1 

8,1 

7,1 

- 6,1 

3,1 

- 4,1 

5 ~f 

2,1 

1,1 


Figure 18. Stiffness Matrix Substitutions for Eight-Node Shear Panel 
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APPENDIX B 

The General Stiffness Equation 


For a shear panel of length A, width B, and thickness T (figure 19a), 
apply a deforming stress a x and a restraining stress o y (figure 19b): 



Figure 19. Unstressed Shear Panel with a x Added 


The unit volume strain energy now stored in the panel is 


u = 





(82a) 


Next, add a deforming stress cr y (figure 20a) 


Ee" 


u = 


X 


2(1 




a y de y 
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Figure 20. Addition of a y and t. Respectively 


Last, by adding a shear stress t (figure 20b)* 


u = 0 ,, E — 57 (e 3 + 2(_ie e + e 2 ) + / rd 7 

2(1 - |i ) x K x y y J 

o 


7 

u ' 2(1 - 7 ) (s x + 2|je x s y + e P + G J 767 


U = 


2(1 - y*) 


+ 2 ^x £ y + S + 2 ^ 2(1 " & 


(82c) 
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Total strain energy is 


A B T 



dx 


2 

e 

x 


+ 2 ^x 





dy dx. 


( 83 ) 


Finally, by using Cas tigliano' s first theorem and taking partials under 
the integral s igns , 


A B 


S = P = 


dU _ ET 


ij i B5 


1 - 1 ? 



Be.. 


x 


e * sT + 11 1 ^ ^ + 6 


'x B5 


Be. 


x 


y B$, 


o o 


+ F ^-2 + I , &L- /, . ,A 

+ e y ^ + 2 y <l •*> 


dy dx. 


(«) 
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APPENDIX C 
Matrix Algebra 


The stiffness matrix [S] of a composite structure, which is also 
positive definite, is inverted to find the^ flexibility matrix [F] by 
first defining a lower triangular matrix [A], i.e., i g j, and its 
transpose such that 


[A] x IA T ] = [S3 

and 

[F] = [A 1 ]" 1 x [A]- 1 , 


or in algebraic notation 



n 

F xj X ^ 

K=i 


(84) 

(85a) 

(86a) 

(86b) 

(87a) 

(87b) 

(85b) 
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For the given matrix [S] of order r, indices are such that X g m g i, 3, 
K, H §■ n ^ r, with those indices representing foundation nodes omitted 
If m = 1 and n = r, then the entire matrix is inverted, otherwise an 
arbitrary submatrix. By definition, an index of zero denotes an element 
of zero value 5 and K f &. 

If nonzero deflections K are desired m addition to or instead of 
a set of loads, then the nodes K are treated as "foundations," and all 
values subtracted from each P 15 such that 


<Vnew = ^old 


S iK 8 K 3 


( 88 ) 


for 1 not a foundation node. The provision is especially useful in 
calculating the stiffness matrix of a complex composite structure of 
shear panels, beams, etc., by systematically assigning all loads zero, 
an arbitrary node a unit deflection, all other nodes zero deflections, 
and applying equations (88), (89), and (90) until the matrix is com- 
plete. Otherwise, once [S] is inverted and {Pj is biased for "nonzero 
foundations," use 


{&) = [F] X {P} 

and 

{P} = [S] x {5} 


(89) 


(90) 


m that order. 
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